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Abstract 

The problem of unicellular-multicellular transition is one of the main issues that is discussing in 
evolutionary biology. In [1] the fitness of a colony of cells is considered in terms of its two basic 
components, viability and fecundity. Intrinsic trade-off function of each cell defines a type of 
cell. We elaborate models providing in [1]. Assuming that all intrinsic trade-off functions are 
linear, we construct a model with different cell types and show that the differentiation of these 
types tends to full specialization. In addition, we attempt to consider the fact that environmental 
factors influence on the fitness of the colony. Thus, we introduce an energy restriction to the 
model and show that in optimum we get situations in which there exists a set of states, each of 
them allowing colony to achieve the same maximum level of fitness. In some states arbitrary 
chosen cell may be specialized, in some - unspecialized, but fecundity and viability of each cell 
belong to limited ranges (which are unique for each cell). It is worth pointing out that the models 
from [1] are not robust. We try to overcome this disadvantage. 

Keywords: unicellular-multicellular transition, differentiation of types, energy constraints, germ- 
soma specialization, life-history evolution. 



1. Introduction 

The problem of unicellular-multicellular transition is one of the main issues that is discussing in 
evolutionary biology. It is necessary to know how colonial organisms are transformed into 
multicellular organisms and what preconditions underlie this transition. 

The separation of a body's tissues is the main characteristic of a multicellular organism. This 
separation means that the majority of cells in the organism is specialized for one function and 
loses the potential ability to specialize for other functions. 

Some cells in a colonial organism may be specialized for specific functions but may not lose the 
ability to specialize elsewhere. If conditions that lead the colony to full specialization are 
performed over a long period of time, it is possible that the unicellular-multicellular transition 
will occur. Therefore, it is important to determine the conditions that contribute to the full 
specialization of a colonial organism. For example, it is reasonable to suspect that different cell 
types may cause full specialization. Another interesting question is the influence of 
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environmental factors on the colony's behavior. We introduce these environmental factors into 
the model in the form of the total energy constraints consumed by considering colony of cells. 

Fundamental models investigated the problem of unicellular - multicellular transition were 
provided in [1]. These models best illustrate evolution of Volvocalean Green Algae but they may 
be applied to other lineages as well. In [1] the model has been illustrated on the example of 
volvocales green algae. These are flagellated photosynthetic organisms with coherent 
glycoprotein cell walls and represent the most appropriate system in research of the process of 
transition under study, because Volvocales linage ranges from the single-cell organisms to 
undifferentiated, soma-differentiated and germ-soma differentiated organisms arranged 
according to the size of the colony. Volvocales live in standing waters and so need flagellar 
beating in order to move toward light and nutrients. Therefore motility is an important factor 
contributing to viability of Volvocales [2]. Volvocales' type of cell division represents palintomy 
with multiple fission. Also, useful fact is that the species with increased cell specialization do not 
have a single origin. 

However, the results in [1] look non-robust. Because of the identity of all cells in the colony 
assumed in [1], in optimum it does not matter which cells belong to the sets of soma- specialized 
cells or germ- specialized cells. So, if we change slightly some characteristics that not reflected 
explicitly in the model, sets of germ and soma-specialized cells changes - the model only 
requires that the ratio between their cardinality should remain constant. Thus, small changes in 
parameters may force soma-specialized cell become germ- specialized immediately. In linear 
case this non-robustness also lies in a fact that no more than the half cells can be soma- 
specialized and no more than the half cells can be germ- specialized. These facts have attracted 
our attention and we provide a new model to overcome above mentioned non-robustness of this 
beautiful model developed in [1]. 

Structure of the text. 

We begin with the model, provided in [1] and describe it in Section 2. Then in Section 3 we 
propose the model with different cell types, extend this model taking into account energy 
constraints in Section 4. In section 5 we give a short survey of related works. Section 6 
concludes. 



2. An overview of the models 

In [1], the authors constructed models that explore the fitness trade-offs at both the cell and 
group levels during the unicellular-multicellular transition. Thus, fitness is considered in terms of 
its two basic components: viability and fecundity. In [1], the trade-off function (1) is studied, 
which reflects the intrinsic relationships that link viability and fecundity within the cell due to 
cell physiology and other constraints. 

Let v be viability and b represents fecundity. Then: 

v = v max -a*b. (1) 
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In [1], the authors noted that in unicellular organisms, the cell must contribute to both of the 
fitness components. In multicellular groups, each cell may be unspecialized, such as in 
unicellular organisms, or, in contrast, may specialize only in the germ or only in the soma. This 
fact can lead to the formation of germ - soma ("G-S") specialization, in which some cells lose 
their autonomy in favor of the group and, as a result, their fitness and individuality are 
transferred from the cell level to the group level. In [1], the cases in which "G-S" specialization 
may occur are studied. It is noted as well that the models are most applicable to volvocine green 
algae. 

There are two types of models that are considered in [1]: the fitness isocline model and the full 
optimization model. We discuss only the second model because this model is more general than 
the fitness isocline model. In the full optimization model, all of the cells are considered 
simultaneously, and the strategic purpose of the colony is to maximize its fitness. Below, we 
describe this model in detail to emphasize its advantages and disadvantages. Then we try to 
improve it. 

2.1. Full optimization model 

Consider a colony consisting of N cells, i = 1 ...N are indices of cells in the colony, b t - 
resulting contribution of cell i to the fecundity of the group, v t - viability-enhancing capability 
of cell i. 

The fitness trade-off function (continuous and determined on a convex hull) is more common 
than a linear function and should be the same for all of the cells from the colony: 

v t = v{b t ). (2) 

The group's level of fecundity is an additive function of variable b;. The group's level of 
viability is an additive function of variable v i5 i.e., 



N N 

B 

i=l i=l 



= Yu bi) v= Yu vi - (3) 



In [1], it is assumed that the group fitness, which we should maximize, is the product of the 
group viability and fecundity: 

W = V * B. (4) 

In [1], the problem of choosing the correct form of the group fitness function is also discussed. 
This type of function (4) is based on simple intuition. For instance, imagine that one cell in the 
group has a high level of fecundity but low viability, and another cell is strictly the opposite, 
with a high level of viability and low fecundity. Each of these cells by itself would have a low 
fitness, but together they can achieve a high fitness for the group. Function (4) considers this 
reasoning in contrast to, for example, the average fitness of all cells (5): 



N 

w= jiY l bi * vt (5) 

i=l 
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However, we should note that most of the assumptions in [1] would still hold even if the fitness 
is described by a more general function with special properties. 



In general, the full optimization model can be written formally as a type of optimization problem 
(6): 



, N N 

W = ^ b t * ^ v t -> max b v 



1 = 1 1 = 1 

Vi = l..N:v t = v(b t ), w 

Assume that there is no initial cost of reproduction in the model. In this case, the following 
results are true: 

1. If the function vQj) is strictly concave, then the group of cells should remain unspecialized. 

2. If the function v(Joi) is linear (v = v max — a * b), then the group of cells behaves as if there 
was just one cell; therefore, each cell is indifferent to specialization. 

3. If the function v(b) is strictly convex, then the group of cells aims for full specialization. In 
addition, if there is an even number of cells in the group, then half should specialize in the germ 

and a half in the soma. If there is an odd number of cells in the group, f^j of these cells should 

specialize in the germ, [^] should specialize in the soma and one cell should remain 
unspecialized. 

An initial investment is necessary for reproduction. This investment requires an additional 
spending of energy, which can lead to the appearance of initial costs of reproduction that can be 
considered within the trade-off function. The initial costs of reproduction lead to full 
specialization in linear and convex cases of improved model and provide the opportunity for 
specialization in the concave case of this model. 

Despite all of the fundamental results of the full optimization model, there are a wide variety of 
problems that cannot be solved using this model and some disadvantages that are connected 
strictly with biological processes that this model cannot describe. For instance, 

1. The full optimization model does not explain why the cells are jointed into a group and exist 
in this state instead of continuing to exist separately. 

2. In the full optimization model the separate form of fitness function of the colony is used. 
However, in some cases, this form is not appropriate. 

3. Models, provided in [1] assumed that all of the cells in the colony should be identical. In some 
cases it is reasonable to omit this identity assumption. 

4. The full optimization model does not strictly reflect the fact that some environmental factors 
may influence on the fitness of the colony. . 



Below we try to construct a model that is based on two important assumptions: 
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1. Cells within the colony are of different types (different intrinsic trade-off functions). 

2. There are some energy restrictions in the system under consideration. 

In Section 2 the linear model, which satisfies only the first assumption, is studied. In Section 3, 
we construct the model that satisfies both of the specified assumptions. 



3. Full optimization model with different cell types 
3.1. Formulation of the problem 

Now we will assume that a set of all cells in the colony is divided into finite number of non- 
intersected, non-empty subsets (we get a partition in the set of all cells), such that: 

- Each two cells belonging to the same subset have the same intrinsic trade-off function (the 
same type); 

- Each two cells from different subsets have different intrinsic trade-off functions (different 
types). 

In order to simplify our analysis we suggest that all of the cells in the colony are unique, so each 
cell has its unique type. Note, that previous, general type of the model can be easily reduced to 
the model with unique cells in linear case. 

Consider a group of N cells, i = 1 ... N - indices of cells in the colony, b t - level of fecundity of 
cell i, Vi — viability of cell i. 

The group's level of fecundity is an additive function of variable b;. The group's level of 
viability is an additive function of variable V;: 

N N 




t=i t=i 



We agree with [1] and continue to apply the type (8) of the group fitness functions because this 
function reflects the synergetic effects of jointly existing cells in the colony, 

W = V * B. (8) 

The main extension of the model from [1] is that each cell has its own parameters of the trade-off 
function. Additionally, we assume that each individual's trade-off function is linear (9). 

Thus, 

Vi = 1JV: Vi = v™ ax - Ui * b u (9) 
Vi = 1JV: 0Ci > 0, 
Vi = 1JV: v™ ax > 0, 
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The main assumption of the full optimization model with different types of cells (strong 
differentiation of types) has following form: 

at ocj, for any i,j G {1, . . , N}, i j. (10) 

We use formula (10) because it presents the pure form of the differentiation of types in the model 
and because models with other types of differentiation can be easily transformed to this model. 

Also, note that everywhere in this work, assumptions (9)-(10) are hold. 

Therefore, we can construct the model as a classic optimization model with constraints, 



= ^ b t * ^ v t -» max bv 



f N N 

W 

\=\ i=l 

i = V j -l "i 

Vi = l,N:bi > 0, 



Vi = 1, JV: Vi = v™ ax - ai * b b W 



'i 

Vi = 1JV: v, > 0. 



This optimization problem can be transformed into a more convenient form (12): 

( N N 

W = y bi* y Vi -» max bv / n n 

^ ^ W = *^« a * -a t *bi) ^max b 



i=l i=l 



Vi = l,N:v i ^v? iax -a i *b i ~ < \ i=1 ;=1 
Vi = 1,JV:6 £ > 0, 



Vi = l,iV:0 <b t <bf ax . 



Vi = 1JV: 17# > 0. 



i 

i=l i=l i=l 

Vi = 1JV: 0 < fc f < b\ 



JV JV JV 

^ bi * (^ i^ ma * - ^ a f * Z7 f ) -» max b 



(12) 



max 



Thus, we should maximize the polynomial of degree two, determined in the hyper parallelepiped 
inR n . 

3.2. An analysis of the model 

In this subsection, we provide some propositions that show the nature and important properties 
of the solution to the optimization problem (12). 

Definition 1. Let T be the class of optimization problems, p G T if and only if p is an 
optimization problem with the form represented below: 



P. n ^ h _ umax 



(13) 



Vi G /£:0 <bi< b\ 



7? c{l,..,7V};/ 2 p c{l,..,tf}; / 3 P ^{l,-.,iV}, 
/[ u/ 2 p u/ 3 p = {l,..,iV}, 
where /[n/ p =/[n/ p =/ p n/ p = 0, (14) 

I 1^1- 

Therefore, any task p from class T has a form (13) and is determined by a triplet of sets 
{jiJ^.'^) satisfying to the conditions (14). For example, problem (12), which is the formal 
description of the full optimization problem with different types of cells, also belongs to class T 
and is matched to the triplet (0, 0, {1, . . , JV}) (denoted as PP). 

Additionally, it is necessary to define some sets that will be useful below: 

F p — represents the domain of optimization task p,p G T, 

H p = [b E R N : 0<bi< b™ ax , Vi G 7 3 p ; b t = b™ ax , Vi G l\; b t = 0, Vi G l\\ V G P, 

K = {H p ,pE P), 

L = {H p> \lV\ = l}. 

{(W (fo) —> TYICLX) 
b G H PP : 31 G L, b is a conditional stationary point of the problem | ^ G j y 

Theorem 1. Let b* represents the solution to the problem (12). Then, the following statement is 
true: 

(b* G Vert{Hpp)) V {b* G F). (15) 
Where Vert(H PP ) represents a set of vertexes of H PP . 
The proof of Theorem 1 is provided in Appendix A. 
Results of Theorem 1 : 

- If all types in the colony are unique, in the optimum there is no more than the one unspecialized 
cell. Furthermore, there cannot be a situation where all of the cells in a group are soma- 
specialized or germ-specialized. Also, this type of model is almost always robust. 

- Consider more general case of this model, where we have a partition in the set of cells and each 
subset of cells from this partition has its unique type (suppose, all subsets from partition have 
cardinalities larger than one). According to Theorem 1, in the optimum some subsets from the 
partition consist from soma- specialized cells; other subsets - from germ- specialized cells; 
finally, the third type of subsets is either empty or contains only one subset of cells indifferent to 
specialization. In biological terms this case describes situations, where some cells in the colony 
are specialized and others established interlayer which is indifferent to specialization and may 
perform different functions depending on different external conditions. This type of model is 
almost always robust too. 
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Furthermore, models provided here allow to more than the half of the cells being specialized on 
soma (or germ) function. 

Thus, we show that cell differentiation in linear cases tends the majority of cells to full 
specialization. 



4. Full optimization model with an energy restriction and different types of cells 
4.1. Formulation of the problem 

All of the previous models have a serious disadvantage: environmental factors do not directly 
influence the fitness of a colony. Suppose that the colony has some way of obtaining energy 
from available resources. Let C be the level of the energy that is available to the colony by using 
a fixed way of obtaining energy and fixed environmental characteristics. It is necessary to have 
ki units of energy in order to construct a unit of fecundity and k 2 units of energy in order to 
support a unit of viability. Considering these judgments, we can construct an energy restriction 
as follows: 

N N 

ki * ^ bj + k 2 * ^ Vi < C, where k x > 0, k 2 > 0, C > 0. (16) 

i=l i=l 

Using this relationship, we obtain a full optimization model with an energy restriction and 
different types of cells. 

( N N 

W = ^ b t * ^ v t -> max b v ' N N 



i=l i=l 

_ „max 



Vi = 1. . N: v t = v™ — cti* bi, 
K * ^ b t + k 2 * ^ vi < C, 

i=l i=l 

Vi = l..N\Vi > 0, 
Vi = l..N:bi > 0. 



W(b) = ^ ^ * ^Jy™* - ai * bi) -> max b 

i=l i=l 

» ^ (17) 

k 1 *2_ i b i + k 2* 2JyT x - cci * bd < C; 

i=l i=l 

Vi = 1..JV : 0 < bi < b™ ax . 



The energy restriction can significantly influence on the behavior of the colony in the model. 
However, assume that this restriction allows colony to exist. 

4.2. An analysis of the model 

Below we provide some propositions that show the nature and important properties of the 
solution to the problem (17). The problem (17) is a task in which we should maximize a 
continuously differentiable function determined on a truncated hyper parallelepiped. 

First of all, it is necessary to introduce sets that will be useful below. 

= C\, 



(N N 
b G R N : k t * ^ ^ + k 2 * Y_t v ? aX ~ a i * b d 
i=l i=l 



8 



IN N 
b G R N : k t * ^ bi + k 2 * ^(v™* ~ a i * b ^ ~ C \> 
i=l i=l 




R = <beH 



Theorem 2. 

Let b* G R N represents the solution to the problem (17). Then, the following statement is true: 

{b* G Vert(H PP ) n M x ) V (b* 6(?)V {b* G R) V [(b* G Vert{H PP n M)) 

A (7? = 0)] (18) 

Where Vert(H PP ) represents a set of vertexes of Z/ PP , 

VertQipp n M) represents a set of vertexes of H PP n M. 

Generally speaking, this theorem determines the location of the optimal point and claims that 
there can be only three cases, each of them can be realized depending on the values of 
parameters of the model: 

1. All cells in the colony are specialized: some in soma, some in germ; 

2. There is the one unspecialized cell in the colony; 

3. There exist a set of states, each of them allows colony to achieve the same maximum level of 
fitness. In some states the cell may be specialized, in some - unspecialized, but fecundity and 
viability of each cell belong to limited ranges (which are unique for each cell). 

We can point out as well that the result of the full optimization model with different types of 
cells is robust. Indeed, small variation of parameters, which are not reflected explicitly in the 
model, cannot lead to sharp changes in the solution. The solution either does not change (as in 
the cases 1 and 2) or change slightly (as in the case 3 due to the fact that the optimal set of states 
is connected). Also, more than half of the cells may specialize in soma or in germ. 

The proof of Theorem 2 is given in Appendix B. 

According to the Theorem 2, in some cases we can get an energetic estimation of the optimal 
value of fitness function: 

Let b* G R N represents the solution to the problem (17). Suppose R * 0. Then, W{b*) > 



4*k 1 *k 2 ' 

meaning that under some circumstances, if the colony has the available level C of the energy, 

c 2 

this energy can be transformed into no less than units of fitness. 

bJ 4*/c 1 */c 2 
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5. A survey of the literature 

There are a lot of papers dealing with the problem of unicellular-multicellular transition in terms 
of fitness and its two basic components - fecundity and viability. Here we provide short 
descriptions of some articles which characterize precisely main developments in discussed 
issues. 

In [2] the problem of transition from unicellular to multicellular organism is studied using some 
physical assumptions concerning different processes in the organism in terms of physical laws, 
for example, hydrodynamic laws. 

In [3] a general mathematical model about the division of labor is introduced. The main 
assumption of the model is that modules can contribute to two different tasks, which connected 
by a trade-off. It is shown that three factors favor that division of labor - positional effects, 
accelerating performance functions and interaction between modules. 

In [4] the models of social choice have been applied to the problem under study. It has been 
shown that applying axiomatic approach allows constructing different social welfare functions 
describing the types of fitness-ranking on the set of alternatives, representing all states of the 
world that are relevant for the group under study. Authors suggest to apply extensive social 
welfare functions and their axiomatic to describe the fitness-functions of colonies and show that 
this axiomatic does not contradict the fact that the transition from unicellular to multicellular 
organism accompanies to replacing concavity by convexity in trade-offs. 

In [5] a diversity in Volvocalean green algae is investigated in terms of the process of transition 
to multicellularity. Authors show that costs of reproduction plays an important role in the 
evolution of multicellularity in Volvocales. The proposed model allows explaining the GS - 
GS/S - G/S form of process of transition to multicellularity. 



6. Conclusion 

We have explored the fitness trade-offs between the basic components of survival and 
reproduction at both the cell and group levels during the unicellular-multicellular transition. We 
have considered the model that describes a colony with different types of cells and linear 
intrinsic trade-off functions. 

As a result, if all types in the colony are unique, in the optimum there is no more than the one 
unspecialized cell. If we have a partition in the set of cells, and each subset of cells from this 
partition has its unique type, then in the optimum one part of subsets from partition consists from 
soma-specialized cells; second part - from germ- specialized cells; the third part is either empty 
or contains only one subset of cells indifferent to specialization. In biological terms this case 
describes situations, where some cells in the colony are specialized and others established 
interlayer which is indifferent to specialization and may perform different functions depending 
on different external conditions. Also, provided models are almost always robust. All in all, we 
show that cell differentiation in linear cases tends the majority of cells to full specialization, 
while in linear cases in [1] the cells are indifferent to specialization. 



10 



Additionally, we add environmental factors which directly influence the fitness of a colony. So, 
we have introduced a power restriction to the model and explored the full optimization model 
with a power restriction and different types of cells. As a result, in optimum we get situations in 
which there exists a set of states, each of them allows colony to achieve the same maximum level 
of fitness. In some states arbitrary chosen cell may be specialized, in some - unspecialized, but 
fecundity and viability of each cell belong to limited ranges (which are unique for each cell). 

We can point out as well that the solution to the full optimization model with different types of 
cells and energy constraint is robust. 



Appendix A. 

Introduce following notation: 

Int A - an interior of a set A, 

Relnt A - a relative interior of a set A, 

OA - a boundary of a set A, 

Vert A - a set of all vertexes of a convex set A, 

V/ - gradient of a function f. 

Proof of Theorem 1. First, let us proof the following 

Lemma 1. For each p G T, | > 1 , let b p G IntF p ; then, b p cannot be a solution to 
optimization problem p. 

Proof. Choose any p G P, > 1. Let b p G IntF p . Suppose b p is a solution to optimization 
problem p. 

1. W p {bp) > W p (b), Vb G Fp, because b p is a solution to optimization problem p. 

2. 3U £ (bp) c F p \ Wp(b p ) > W p (b),Vb G U £ (b p ), because b p is a solution to optimization 
problems p and b p G IntFp. Therefore, b p is a local extremum point of function W p . 

3. Because W p is a continuously differentiable function on the entire domain, the necessary 
condition for a local extremum of the function W p at b p has the form: VW p {bp) = 0. 



VW p (bp) = 0 

dW 



d -^{b p ) = [y v r x + £ »r x aj * ( £ br x + £ t t ) = o,v ; - g / 3 p ; 

; Kiel? ie/ 3 p ie/ 3 p / \ie/f ie/ 3 p / 
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^ (a t + a,-) * bf = I ^ i^"* + ^ ~ a i * Z ^ ) ' V; G ^ 

iEl v 3 \ie/ 2 p ie/ p ie/ P / 

We obtain the system of | linear equations. Since |/g | > 1, we can choose two different linear 
equations with indexes j and k accordingly: 

£(a, + a y ) * ftf = ( £ ^ + £ vr x - «/ * £ *™ 

ie/ p \ie/ 2 p te/ p ie/ P 

^ (a, + a k ) *bf = 1^ v ? aX + X ^ ~ ak * Z ^ ) 
.ie/ p \ie/ p te/ p te/ P / 

Subtract the second equation from the first equation. Since 

£ (a ; - a k ) * bf = -(ay - a fc ) * £ 



then 



te/ 3 p te/£ 



\ie/ p iE/ P / 
^ fcP = - ^ < 0 -> 3i G 7 3 p : fcf < 0 -> 6 p g F p 



te/ P te/ P 



We have a contradiction; therefore, b p cannot be a solution to optimization problem p. Lemma is 
proved. 

Now we prove Theorem 1. Consider the optimization problem (12). 

1. If N=l, then, obviously, the statement (15) is true. 

2. Suppose N>1. Then, optimization problem (12) belongs to the class T. Consider a point b G 
Int(Hpp), (note, that H PP = F PP ). According to Lemma 1, b cannot be a solution to optimization 
problem (12); therefore, the solution to (12) belongs to the set dH PP . 

Consider the behavior of function W in 9H PP . Choose one variable and prescribe it the value 0 
or bf 135 '. Thus, we describe a set dH PP = U pe j>.|i p |=n-i ^p- Therefore, we obtain tasks, each of 

which would be equivalent to some pG Pil^l = N-1. According to Lemma 1, for each 
task p G J 3 : 1 13 I = N — 1, for all b G ReInt(H P ), it follows that b is not a solution to the problem 
(12). Therefore, we should find the solution to (12) in the set U pEy .| I P| =N _ 2 H P where our 

conclusions are similar. We should repeat our procedure until we obtain tasks in the edges and 
vertexes of H PP . The solution to (12) belongs to the vertexes or edges of hyper parallelepiped 
(12). Therefore: 

12 



(b* G Vert(Hpp)) V (b* G F), 

So, Theorem 1 is proved. 



Appendix B. 

Proof of Theorem 2. First, we state the lemma. 
Lemma 2. Consider optimization problem (19): 

N N 



i=l i=l 
N N 



(19) 



k x * ^ b t + k 2 * Yf?™* ~ «i * b t ) = C, 

i=l i=l 

Vi = l..N:0<b t < b™ ax . 

Let b** G R N be the solution to problem (19). Then, the following statement is true: 

[b** G R] V [(V* G Vert{H PP n M)) A (R = 0)] (20) 

Proof. First, look at the class of optimization problems 11 , which is constructed based on the 
coefficients of the task (19). Each task u(I?,I$,If) Ell has a type (21) (where I? Q 
{l,..,Nl q c {1, . . , N], If c {1, ..,7V} : 7 P U 7 P U 7 P = {1, . . , JV}; /? n 7 P = 7 P n /f = 7 P n 
7 P = 0; |7 3 w | > 2): 



ki * ( Z ^ + Z ^ ) + kz * ( Z ^ + Z ^ " Z ai * bi ) = c (2l-2) 



Each u(7j\7", 7") G 17 represents a task of maximizing the function of the appropriate form 
(21.1) on the hyper plane of the appropriate form (21.2) in j?' 7 "' space. We choose and fix some 

u(7", 7", 7") = uQ lt 1 2 , h) G ^ an( l simplify the form of optimization problem u. 

(B(bi\i G 7 3 ) * V(bi\i G 7 3 ) -> max 6 . iEf3 (22ll) r22^ 
Ifci *B{b t \i G 7 3 ) +/c 2 *V r (fc i |i G 7 3 ) = C (22.2) ^ J 

1 

w(B(6 £ |i G 7 3 )) = — * B(bi\i G 7 3 ) * (C - k t * S(7) £ |i G 7 3 )) -» max fc(>i6 , 3 
fe 1 *B(fc t |ie/3) + fc 2 *7(^|i6/ 3 ) = C 
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B{b i \iei 3 ) = \ybr x +Y J b i 



where: < 



te/i iel 3 



vm £ / 3 ) = I £ + X - Z ai * ^ 

\ie/ 2 ie/ 3 ie/ 3 , 



This task is related to the problem below: 
f(bi\i G 7 3 ) = w(B(bi\i G 7 3 )) = ^B(bi\i G 7 3 ) * (C - k t * B{b t \i G 7 3 )) max 6 .; e/3 (23) 
The necessary condition for a local extremum of the function / has the form: 
dw\B(b 1 ;...;b\ h \)) _ 9w ( 5 ) dB(b t \i G 7 3 ) _ dw(S) 

This system is equivalent to one equation: 



= 0; V i G 7 3 



dw(5) C 

Moreover, it is important to note that the function f(bi\i G 7 3 ):7?' /3 ' -> 7? is concave inT?'' 3 ', 
because: 

1. The function w(5) = — 5 * (C — /q * 5) is concave in R. 

2. The function B^; 6|, 3 |) = (I ie/i frf 1 "* + Ete/ 3 ^t) is linear in its domain. 
Suppose x G 7? |/3, ,y G 7? |/sl . Then: 

f{X * x + (1 - X) * y) = w(5(l * x + (1 - X) * y)) = w{X * 5(x) + (1 - X) * 5(y)) 

> X * w(5(x)) + (1 - X) * w(5(y)) = X * f{x) + (1 - X) * f{y), Vl G [0; 1] 

These equalities indicate that each stationary point of the function/"^; ...;b\ l3 \) in7?' /3 ' space 
represents the point of a local and a global maximum of the function fib-^; ... ; fr|/ 3 |). Moreover, 
all of the points of the global maximum of / belong to a (|7 3 |-l)-dimensional hyper plane, which 
is described by the equation below: 



2> = 27^ -Z 6 '"*" < 25 > 

t6/, iei 



Task (22) is equivalent to task (23) with restriction (22.2). This restriction represents a (|7 3 | — 
1) -dimensional hyper plane as well. As a result, all of the points that belong to the intersection of 
planes (22.2) and (25) are solutions to a problem of maximizing function (22.1) with restriction 
(22.2). 

Therefore, the solutions to problem (22) satisfy the system of two linear equations: 
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' * ymax 



- k 2 * adb t = C - h * £ -*2 * ( £ *r" + ^ v, 

V.tE/ 3 iE/i \te/ 2 i6/ 3 



(26) 



max 
i 








J = a* 









These two planes are parallel or represent one plane when their normal vectors are linear- 
dependent vectors: 



-» k 1 -k 2 *a 1 = k ± -k 2 * a lhl -> a x = a lh 



According to assumption (10) about values of parameters of the model, a x oc\i 3 \. As a result, 
planes under consideration have a nonempty intersection; moreover, solutions to problem (22) 
belong to a (|/ 3 | — 2) -dimensional plane in i?'' 3 ' space that can be described using system (26). 
Additionally, it is important to mention the statement below: 

Let b satisfy (26). Then: 

WQ>) > W(b x ), Vb x G M, but not satisfying (26). (27) 

This proposition follows from the curvature of optimized function and its concavity. 

Note that all of the statements considered earlier in this lemma for the fixed optimization 
problem uQ lr I 2 , 1 3) remain valid for each u(7", - ) £ V~ 

Now we have a useful approach to investigate the problem (19). We emphasize the fact that the 
problem (19) represents the task u(0, 0, {1, . . , iV}) G U with additional restrictions 0 < b t < 
ffnax yj _ i n j± s f oun( i earlier, all of the solutions to task u belong to a (N-2)-dimensional 
plane, that is why there are two possible cases: 

1. There exists at least one point that simultaneously belongs to the optimal hyper plane and to 
the hyper parallelepiped H PP = {b G R N : 0 < b t < bf ax ,Vi = 1..N}. According to condition 
(27), these and only these points represent the solutions to (19). 

2. There is no point that simultaneously belongs to the optimal hyper plane (26) and to the hyper 
parallelepiped H PP . According to condition (27) and due to the concavity of optimized function 
in energy restriction, we should search for the solution to the task (19) in dH PP n M. We also 
obtain tasks u(/", /",/") G 11, which are determined in appropriate hyper parallelepipeds. 
Because initial system (26) has no solution, corresponding systems (26) constructed for reduced 
tasks would have no solutions. We can repeat this procedure until we receive vertices of the 
set H PP n M. 

Therefore, we show that the following statement is true: 

[b** G R] V [0** G Vert{H PP n M)) A (i? = 0)]. 
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So, Lemma 2 is proved. According to Theorem 1 and Lemma 2 it is obvious that the following 
statement is true: 

{b* G Vert(H PP ) n M x ) V (6* 6(?)V [6* G R] V [(V G Vert(H PP n M)) A (7? = 0)]. 
Theorem 2 is proved. 
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